From scattering theory to complex wave dynamics in non-hermitian PT-symmetric resonators 
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I review how methods from mesoscopic physics can be applied to describe the multiple wave scattering and 
complex wave dynamics in non-hermitian PT-symmetric resonators, where an absorbing region is coupled 
symmetrically to an amplifying region. Scattering theory serves as a convenient tool to classify the symmetries 
beyond the single-channel case and leads to effective descriptions which can be formulated in the energy domain 
(via Hamiltonians) and in the time domain (via time evolution operators). These models can then be used 
to identify the mesoscopic time and energy scales which govern the spectral transition from real to complex 
eigenvalues. The possible presence of magneto-optical effects (a finite vector potential) in multichannel systems 
leads to a variant (termed VTT' symmetry) which imposes the same spectral constraints as VT symmetry. I 
also provide multichannel versions of generalized flux-conservation laws. 



I. INTRODUCTION AND MOTIVATION 

Optical systems combining lossy and active elements with 
carefully balanced absorption and amplification rates provide 
a platform to implement analogues of non-hermitian VT- 
symmetric quantum systems [[IKS. This allows to realize 
unique optical switching effects Q-y/ZD which are based on 
the transition in the spectrum from real to pairwise complex- 
conjugate eigenvalues. 

The wave scattering and dynamics in these systems provide 
an aspect which also permeates mesoscopic phenomena in mi- 
crocavity lasers II 1811. coherent electronic transport 11191 12011. 
superconductivitv |[l9l l2lh . and quantum-chaotic dynamics 
Il22ll . The main goal here is to relate how tools established 
in these disciplines can be used to approach the exciting spec- 
tral and dynamical features in PT-symmetric optics. Scatter- 
ing theory [0j|[2(l[23|], in particular, is ideally suited to deal 
with the complications of complex wave dynamics in multi- 
channel situations, which only appear when one goes beyond 
one-dimensional situations, and also easily accounts for any 
leakage if a system is geometrically open, as is often required 
by the nature of the desired optical effects, or because of the 
design of practical devices. 

We thus describe in detail how this approach can be applied 
to optical realizations of PT-symmetric systems with a wide 
range of geometries lfl3l - [l7[ 0-[26|]. (For PT-symmetric 
scattering in effectively one dimension see, e.g., Refs. 11271- 
l29h .) To set up the approach under these conditions it matters 
whether a geometric symmetry inverts or preserves the hand- 
edness of the coordinate system, especially in the presence 
of magneto-optical effects (vector potentials, which generally 
cannot be gauged away beyond one dimension). An analo- 
gous bifurcation arises in the specification of the time-reversal 
operation. Thus the formulation of VT symmetry itself re- 
quires some care. In particular, an alternative appears, termed 
VTT' symmetry II24I1 . which results in the same spectral con- 
straints but imposes different physical symmetry conditions. 
The scattering approach also allows to restate the symmetry 
requirements in terms of generalized conservation laws Il6h . 
which we here extend to multiple channels. Furthermore, the 
approach leads to effective models of complex wave dynam- 
ics which relate the spectral features to universal mesoscopic 
time and energy scales l24l - [26ll . 



We start this exposition with a brief recapitulation of the 
analogy between optics and non-hermitian quantum mechan- 
ics (section |II]), and discuss variants of geometric and time- 
reversal operations which can be used to set up VT and VTT' 
symmetry (section UTT] and section HV1). We then describe how 
scattering theory can be applied to study the spectral features 
in these settings, first generally (section |V]-section IVIII) and 
then for a coupled-resonator geometry (section I Villi) . This 
leads to effective models (sectionllXt which can be formulated 
in the energy domain (via Hamiltonians) and in the time do- 
main (via time evolution operators), and capture the relevant 
energy and time scales (section[X]). The concluding section IXj 
also describes possible generalizations of these models. 



II. OPTICAL ANALOGUES OF NON-HERMITIAN 
QUANTUM MECHANICS 

The PT-symmetric optical systems mentioned in the in- 
troduction cover a large range of designs, including coupled 
optical fibres, photonic crystals and coupled resonators, as de- 
picted in figure [T] Some of the elements are absorbing while 
others are amplifying, and the absorption and amplification 
rates, geometry, and other material properties are carefully 
matched to result in a symmetric arrangement. 

We concentrate on two common designs, effectively two- 
dimensional systems (planar resonators or photonic crystals, 
with coordinates x, y) as well as three-dimensional pillar-like 
systems (such as arrangements of aligned optical fibers with 
a fixed cross section) where the geometry does not depend on 
the third coordinate z, and assume that the relevant effects of 
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FIG. 1: (Online version in colour.) Common designs of PT- 
symmetric optical systems with absorbing (light green) and amplify- 
ing (dark red) elements, (a) Two optical fibers, (b) Localized modes 
in a photonic crystals, (c) Coupled-resonator geometry. 
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wave scattering, gain and loss can be subsumed in a refrac- 
tive index n(r), where Im n(r) > signifies absorption (loss) 
while Imn(r) < signifies amplification (gain). Later on, 
we also include a vector potential A(r) representing magneto- 
optical effects. 

These assumptions allow to separate transverse magnetic 
(TM) and electric (TE) modes, with the magnetic or electric 
field, respectively, confined to the xy-plane (thus transverse to 
z). The electromagnetic field can then be described by a scalar 
wave function ip(r), which represents the ^-component of the 
electric or magnetic field, respectively. One now arrives at two 
principal situations, which both serve as an optical analogue 
of non-hermitian quantum mechanics. 



A. Helmholtz equation 

On the one hand, one can focus on the propagation in the 
xy-plane, resulting in an effectively two-dimensional system 
which is described by a Helmholtz equation 




£(uj) = 



v = 



(1) 



The Helmholtz equation is analogous to a stationary 
Schrodinger equation, thus, an eigenvalue problem for the 
frequencies uj, which can become complex due to the loss 
and gain encoded in n, or because of leakage at the bound- 
aries of the system in the propagation plane. This then de- 
scribes quasistationary states which decay (Im uj < 0) or grow 
(Im uj > 0) over time. 



B. Paraxial equation 

On the other hand, one can focus on the propagation into 
the z-direction perpendicular to this plane, which is then often 
described by a paraxial equation 




£(«Mr) = 0, 

v 2 + " 2w 2 2 w 



^(r)V^V 



fir 

u; 2 n 2 (r) 



(TM) 

(TE) 



(2) 



Here ip(r) now is an envelope wave function, obtained after 
separating out a term exp(iftz). This equation remains valid 
if n varies only slowly with z, \d z n/n\ <C 

The paraxial equation is an analogue of a time-dependent 
Schrodinger equation, thus, a dynamical equation where an 
initial condition is propagated forwards — here, not in time, 
but along the spatial direction with coordinate z. It is then nat- 
ural to probe the system at an initial and a final cross-section, 
while the physical frequency uj is now real. 

As we assume that n is ^-independent, the paraxial equation 
is associated with the eigenvalue problem C{k)^{t) = for 
the propagation constant k. This eigenvalue problem shares 
all relevant mathematical features with the problem (Q]) for 



uj (in particular, the eigenvalues again can be complex be- 
cause of loss, gain, and leakage in the system); so do variants 
that rely on wave localization in the cores of optical fibers or 
Wannier-like modes in a photonic crystal, with the gradient 
replaced by hopping terms. For definiteness, we shall employ 
notations adapted to the eigenvalue problem for uj. 



III. SYMMETRIES OF THE WAVE EQUATION 

The investigation of symmetries in higher-dimensional sys- 
tems of arbitrary geometry is a central aspect of mesoscopic 
physics. For hermitian systems, a complete classification re- 
quires to take care of magnetic fields, internal degrees of free- 
dom such as spin, pairing potentials and particle-hole symme- 
tries fl9U22ll . Non-hermitian systems provide an even wider 
setting, with a mathematical classification, e.g., provided in 
Ref. 1 30]. In order to identify and specify the role of sym- 
metries for their optical realizations we denote the operator in 
the wave equation by C(uj; n(r)), which explicitly takes care 
of the functional dependence on the refractive index n. 



A. VT symmetry 

In VT- symmetric quantum mechanics fl-d], the parity op- 
erator V generally stands for a unitary transformation which 
squares to the identity. In optical systems, this operation is 
usually realized geometrically via an isometric involution, still 
denoted as V, which inverts one, two, or three coordinates. To 
specify the consequences for the wave equation we promote 
V to a superoperator (we employ this notation as we will also 
encounter a transformation T' which does not correspond to 
an ordinary unitary or antiunitary operator). Then 



V[C(uj ] n(r))]=C(uj;n(Vr)), 



(3) 



which should be read as a rule how to write the wave equa- 
tion in the transformed coordinate system. The transformed 
equation is then solved by Vip(r) = ip(Vr). 

Conventional time reversal is implemented by complex 
conjugation in the position representation, which constitutes 
an antiunitary operation. We then have 



T[C(uj;n(r))]=£(uj*;n*(r)), 



(4) 



which delivers a wave equation solved by Tip(r) = ^*(r). 

The wave equation now displays VT symmetry if the re- 
fractive index obeys n(r) = n*(Vr). In this situation V is an 
involution that interchanges amplifying and absorbing regions 
with matching amplification and absorption rates, and 

PT[£( W ; n(r))] = C(w* ; n* (Vr)) = C(co*; n(r)), (5) 

which constraints the spectral properties of C if the boundary 
conditions also respect the symmetry (this is typically not the 
case if the system is open). The eigenvalues uj n then obey 



UJ n LJ= 



(6) 



and thus either real (n = n) if -0* (Vr) linearly depends on 
ip n (r), or occur in complex-conjugate pairs (n ^ n) if that is 
not the case. 

B. T symmetry 

It is important to distinguish the conventional antiunitary 
time-reversal operation T from another operation that is also 
often termed time-reversal (HQ. If £(r) were hermitian (uj 
and n both real with appropriate boundary conditions), then 
the action of T could not be distinguished from the action of 
a superoperator T' acting in the position representation as 

T'[C(uj;n(r))} =£ T ( W ;n(r)). (7) 

Thus, T' transforms the right eigenvalue problem Cip = 
into the left eigenvalue problem C T ip = 0. This delivers the 
same spectrum, and while the right and left eigenfunctions ip 
and ip for a given eigenvalue generally differ they are related 
by biorthogonality constraints. In the presently assumed ab- 
sence of magneto-optical effects, T' is indeed an exact sym- 
metry, C t (uj; n(r)) = C(uj;n(r)) and thus ip = This 
holds even if n is complex, thus, even if T symmetry is bro- 
ken. In the non-hermitian case, therefore, these two operations 
are distinct and must be treated separately. 

IV. MAGNETO-OPTICAL EFFECTS 

In mesoscopic systems with complex wave dynamics, the 
alternative time-reversal operation T' governs a multitude of 
effects ranging from coherent backscattering and wave lo- 
calization to minigaps in mesoscopic superconductors fl9l - 
|2ll|3il]. To clarify the role of this operation we now consider 
magneto-optical effects, which are described by a (possibly 
complex) vector potential A(r) that enters the wave equa- 
tion through terms of the generic form (V + iA(r)) 2 . For 
this purpose we denote the operator in the wave equation as 
£(cj;n(r), A(r)). 

A. VT symmetry 

The involution V transforms the vector potential according 
to (V + zA(r)) 2 -> (VV + iA(Vr)) 2 = (V + i[VA](Vr)) 2 , 
as V is an isometry. Thus 

P[£( W ;n(r),A(r))]=£( W ;n(Pr),[PA](Pr)). (8) 

The resulting magnetic field depends on whether V preserves 
or inverts the handedness of the coordinate system, as V x 
VA(J>v) = det(V)V[V x A(r)]. 

The time-reversal operator T transforms (V + iA(r)) 2 — >> 
(V-zA*(r)) 2 ,so 

T[£(a;;n(r),A(r))] = n*(r), -A*(r)). (9) 

VT symmetry thus requires 

n(r) = n*(Vr), A(r) = -VA*(Vr). (10) 
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FIG. 2: (Online version in colour.) Illustration of magneto-optical 
effects (associated with an effective magnetic field B) for (a) VT 
symmetry and (b) VTT symmetry, where V is a reflection l24f] . 

E.g., when V is a reflection x — >• — x this holds for a homoge- 
neous magnetic field that points in the z direction [see figure 
Ha)]. 

B. VTT symmetry 

We now identify a modified symmetry, which eventually 
yields the same spectral constraints as obtained for VT sym- 
metry, but imposes a different condition on the magnetic field 
i24ll . This variant follows from the inspection of the T' oper- 
ation. In the position representation d J = —d x ,dy = —d y , 
and thus 

T f : (V + iA(r)) 2 [(V + iA(r)) 2 ] T = (V - iA(r)) 2 . 

(11) 

This effectively inverts the vector potential, 

T'[C(lo; n(r), A(r))] = C(co; n(r), -A(r)), (12) 

which confirms that T' symmetry is broken when A is finite. 
Consider now 

VTT'[£(uj; n(r), A(r)] = £(w*;n*(Pr),PA*(Pr)). 

(13) 

This turns into a symmetry if 

n(r) = n*(Vr), A(r) = VA*(Vr). (14) 

E.g. , when V is a reflection x — >> — x this holds for an antisym- 
metric magnetic field pointing in the z direction [see figure 
12b)]. 

V. SCATTERING FORMALISM 

The optical systems considered here are naturally open, 
with leakage occurring, e.g., at the fiber tips and waveguide 
entries, or because the confinement in the cross-section relies 
on partial internal reflection at refractive index steps or semi- 
transparent mirrors. These systems can thus be probed via 
scattering (El H(| 0L which delivers comprehensive insight 
into their spectral and dynamical properties and illuminates 
the consequences of the symmetries discussed above, as well 
as the role of multiple scattering addressed later on. 

These properties are encoded in the scattering matrix S(u), 
which relates the amplitudes ai n?n , a out , n in incoming and 
outgoing scattering states Xm,n, Xout.n, 

a ou t = 5(a;)a in . (15) 
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This relation is generally obtained by the solution of the wave 
equation under appropriate conditions at the boundary dfl 
of the scattering region ft (outside of which we set n = 1, 
A = 0). The scattering states are assumed to be flux orthonor- 
malised, 

/ dS- [x* n Vxa',m-X<r',mVx* n ] = 2iaS nrn S (J(jf , (16) 
Jon 

where a = 1 for outgoing states and a = — 1 for incoming 
states. Two popular choices are states with fixed angular mo- 
mentum in free space, and transversely quantized modes in a 
fixed- width waveguide geometry. 

We now describe the adaptation of this formalism to VT 
and PTT'-symmetric situations jTM^HHH. 



A. VT symmetry 

The symmetries of a scattering problem are exposed when 
one conveniently groups the scattering states. To inspect V we 
call half of the incoming states 'incoming from the left', and 
the other half 'incoming from the right' . This does not need 
to be taken literally; all what matters is that the two groups are 
converted into each other by V. The same can be done for the 
outgoing states. The scattering matrix then decomposes into 
blocks, 

where r describes reflection of left incoming states into left 
outgoing states, r' describes the analogous reflection on the 
right, while t and t' describe transmission from the left to the 
right and vice versa. The parity V interchanges the amplitudes 
of the left and right states, which can be brought about by a 
a x Pauli matrix, 

V[S] = a x Sa x = (J,' (18) 

For consideration of the T operation, we group incoming 
and outgoing states into time-reversed pairs. When T acts on 
a wave function the amplitudes become conjugated, while the 
frequency in the wave equation changes to cj*. Furthermore, 
incoming states are converted into outgoing states, so that the 
relation (TT3T) must be inverted. Thus, 

T[S(lo)} = (19) 

In combination, we have 

VT[S(oj)}=a x {S- 1 (uj*)}*<T x . (20) 

A system with VT symmetry is then characterized by the in- 
variance 

* x {S- 1 (uj*)}** x = S{u>), (21) 



which results in the constraint a x S*(u*)a x S(uj) = 1, or 
t f *(uj*)r(uj)+r*(uj*)t(uj)=0, 

r*(u>*)r'(u))+lf*(u)*)lf(u)) = 1, 
r'*(u*)r(uj) + t*(u*)t(u) = t. (22) 



B. VTT' symmetry 

In hermitian problems the scattering matrix is unitary if uj is 
real, and the T operation is equivalent to the operation i[l9l l20ll 

T'[S(u)}=S t (uj), (23) 

which corresponds to the operation identified above by inspec- 
tion of the wave equation. In non-hermitian settings, this de- 
livers the solution of the scattering problem associated with 
the transposed wave equation C T ip(r) — 0, 

a out =S T (uj)8L in . (24) 

In ordinary optics, where A = 0, T' remains an exact sym- 
metry even when n and uj are complex. We then find the im- 
portant constraint S(uj) = S t (uj), thus 

r = r T , r = r /T , t' = f '. (25) 

The combined operation 

Prr[S(oj)] = a x {S-\u;*)^a x (26) 

turns into a symmetry if i24ll 

t f (cj*)r(cj) + r f (cj*)t(a;) = 0, 

r'\u)*)r(ui) + t'\uj*)t(uj) = 1. (27) 

This is realized for systems obeying the symmetry require- 
ments ([TQb . 



VI. GENERALIZED FLUX CONSERVATION 

For hermitian systems (uj, n and A all real), the unitarity 
S^S = 11 of the scattering matrix ensures the conservation of 
the probability flux, a? n ai n = aQ Ut a out . In Ref. fl6ll analo- 
gous conservation laws were established for one-dimensional 
non-hermitian systems with VT symmetry, and it was found 
that these laws are automatically guaranteed by the symmetry 
conditions (l22l) . These flux conditions only apply at real uj, but 
provide a useful alternative perspective on the role of symme- 
try. Here we extend these conditions to higher-dimensional 
(multichannel) systems, include magneto-optical effects, and 
also allow for VTT' symmetry (these considerations are orig- 
inal). 
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To illustrate the general strategy we consider the Helmholtz 
equation for TM polarization, with A = and uj real. 
Thus ip(r) and ip*(Vr) both solve the same wave equation 
£(cj)^(r) = £(uj)i/j* (Vr) — 0, with £ specified in equation 
©. Now take the following volume integral over the region 

a 

= / dr[^(Vr)£(uj)t/j(r) -^(r)£(uj)^(Vr)} 

= [ dr[^(Vr)V 2 ^(r)-^(r)V 2 ^(Vr)} 

= [ drV[^(Vr)V^(r)-^(r)V^(Vr)} 

= [ dS-[^(Vr)V^(r)-^j(r)V^(Vr)], (28) 
Jon 

under application of Stoke 's theorem. The resulting surface 
integral can be evaluated with equation (QjJ). This delivers the 
generalized flux-conservation law 

a^cr^ain = a^cr^aout = s.}^ (u)a x S(u)R in , (29) 

which together with S(uj) = S t (uj) (as A = 0) amounts to 
the constraints (I22L specialized to the case where uj is real. 

The same constraints follow for TE polarization (as n = 1 
outside Q, so that the refractive index does not feature in the 
surface integral). Next, we include a finite vector potential A. 
For VTT' symmetry, one can follow the steps given above, 
which directly delivers the constraints (f2Tb (again specialized 
to real uj). The case of VT symmetry with finite vector po- 
tential is more intricate, as one then needs to invoke the trans- 
posed wave equation £ T (uj)i()(r) = £ T (uj)ijj* (Vr) = 0. One 
then integrates 

0= f dr[ij*(Vr)£(uj)iP(r) -iP(r)£ T (uj)ij*(Vr)] 

= [ dS-[^*(Pr)V^(r)--0(r)V^*(Pr)], (30) 
Jan 

which results in the generalized flux-conservation law 
a^a^ain = a.l ut a x s. out . In combination with equation (l24k 
this condition now requires cr x = S*(uj)(j x S(uj), which is 
again automatically fulfilled if the scattering matrix obeys the 
symmetry constraints (l22l) . 

In all these cases, the generalized flux-conservation rela- 
tions at real uj thus does not impose any extra conditions be- 
yond the symmetry requirements of the scattering matrix. 

VII. SCATTERING QUANTIZATION CONDITIONS 

It is now interesting to ask how the spectral properties of 
VT or PTT'-symmetric systems emerge within the scatter- 
ing approach. We contrast the case of open systems (where the 
symmetry operation relates quasibound states with perfectly 
absorbed states, thus, connects states of a different physical 
nature) with closed systems (where the symmetry operation 
relates ordinary bound states, thus, connects states of the same 
nature, whose spectral properties then are constrained). 



A. General quantization conditions for quasibound and 
perfectly absorbed states 

Quasibound states fulfill the wave equation with purely out- 
going boundary conditions, thus, ai n = but a out finite. In 
view of equation ([T5l) , this requires 

||5H|| = oo, (31) 

and so uj has to coincide with a pole of the scattering matrix. 
This results in a quantization of the admitted frequencies uj n , 
which in general are complex. 

Quasibound states describe systems whi ch gene rate all ra- 
diation internally, and in particular, lasers Ill3l4l7ll . In a pas- 
sive system the poles are confined to the lower half of the 
complex plane, with the decay enforced by the leakage. In an 
active system, however, gain may compensate the losses and 
result in a stationary state, which signifies lasing. The thresh- 
old is attained when the first pole reaches the real axis. 

Recent works have turned the attention to perfectly ab- 
sorbed states, for which the role of incoming and outgoing 
states is inverted EHED]- These boundary conditions trans- 
late to 115(^)11 = 0, which results in quantized frequencies 
uj n coinciding with the zeros of the scattering matrix. 

For a passive system, where TT'[£(uj)] = £(uj*), a quasi- 
bound state ip n (r) can be converted into a perfectly absorbed 
state i/j n (r) = ijj* (r) by the TT' operation [for exact T' sym- 
metry one simply has ip n (r) = ^* (r)]. This state then fulfills 
the wave equation at uj n = cj*, which confines the zeros to 
the upper half of the complex plane. 

For non-hermitian systems the relation between the poles 
and zeros is in general broken. With VT or VTT' invariance, 
however, the states ^ n (r) = i/j^(Vr) or i/j n (r) = ip*(Vr) are 
paired by the respective symmetry. Moreover, the frequencies 
uj n = cD* then are no longer constraint to the lower half of 
the complex plane. At the lasing condition uj n = 0, the las- 
ing mode is then degenerate with a perfectly absorbed mode, 
which leads to the concept of a PT-symmetric laser- absorber 



B. Closed resonators and spectral constraints 

By taking the appropriate limit of the expressions for open 
systems, the scattering approach to mode quantization can be 
used to study closed systems, for which VT symmetry was 
originally defined. The quasibound states then turn into nor- 
mal bound states, and one recovers the spectral constraints ©. 
These considerations also apply to VTT' symmetry. They 
can also be based on the perfectly absorbed states, which be- 
come degenerate with the quasibound states. 

When one reduces the leakage from a non-hermitian VT 
or VTT' -symmetric optical system, it will ultimately start to 
lase. If some frequencies of the closed system are complex 
then lasing is attained at a finite leakage otherwise 
one approaches at-threshold lasing Il3h . 
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FIG. 3: (Online version in colour.) Scattering amplitudes and their 
relation by reflection and transmission blocks of the scattering ma- 
trix, for the coupled-resonator geometry considered in section IVlTll - 
section[X] In (b), a finite transparency T of the interface is taken into 
account. See also 



VIII. COUPLED-RESONATOR GEOMETRY 

The economy of the scattering approach is full exposed 
when one applies it to multichannel systems capable of dis- 
playing complex wave dynamics. Here we review this for 
a broad class of systems in which an absorbing resonator is 
coupled symmetrically to an amplifying resonator, as shown 
in figuieElEilBlMlll. 



For a PT-symmetric system obeying equation (f22l) this con- 
dition takes the form 

detKH-IrU^)-^^*)^ 1 ^*)^^*)]*) =0, (36) 

while for a PTT-symmetric system obeying equation (f27l) 
this gives 

det(^M-[^(^)-t i (o;*K 1 ( W *)4(a;*)]t) = 0. (37) 

For a closed system the scattering matrices reduce to the 
reflection blocks Sl{u) = r' L (oS) and S 'r(oj) = r r(uj) , while 
the transmission vanishes. In the case of VT symmetry, with 
tr(u) — {[ r L(^*)] -1 }*> equation (l36b assumes the form 



det(riH-[r' L (a;*)]*)=0, 



(38) 



which entails the constraints ©. On the real frequency 
axis, this reduces to the condition det lmr f L (uj) = 0. Analo- 
gous observations also hold for closed VTT'- symmetric res- 
onators, for which the quantization condition (ITTt becomes 



det(ri>)-K( W *)]t)=0. 



(39) 



A. Scattering matrix 



Under inclusion of a semitransparent interface, similar con- 
ditions can be derived from the general expression 



Denoting the scattering matrices of the two resonators as 



S L 



tL r' 



Sr 



r R tR 
t R r R 



(32) 



the scattering matrix S = Sl ° Sr of the composed system is 
given by 




-r R t L 



t' 1 t' 

L L l-r R r' L R 



tRj 



^—r't' 

-r L rR L tt 



. (33) 



VT symmetry follows if the scattering matrices are related by 
Sr = VT[Sl], with this operation specified in equation d20lh 
while VTT' symmetry is realized if the scattering matrices 
are related by Sr = VTT' [Sl], as specified in equation d26l) . 

This construction can be amended to include an interface 
of finite transparency T, with the scattering matrix, e.g., spec- 
ified by 



S T = 



iVr 



(34) 



The total scattering matrix of the system is then given by S = 

Sl o S t ° Sr. 



B. Quantization conditions and spectral constraints 

With equation (l33l) the scattering quantization condition 
(l3Tt becomes 



det [ri(4^)-l]=0. 



(35) 



det 



r' L 
r R 



0. 



(40) 



These scattering quantization conditions can all be inter- 
preted as conditions for constructive interference upon re- 
turn to the interface between the two parts of the resonator. 
The analogous conditions for perfectly absorbed states follow 
from the replacement uj — » uj*. 



IX. EFFECTIVE MODELS 

The coupled-resonator geometry allows to make contact to 
well- studied standard descriptions of multiple scattering lfl9i 
l23ll33l - [35ll . This leads to effective models of complex wave 
dynamics in systems with VT and VTT' symmetry, which 
we first formulate in the energy domain [|24j1. and then in the 
time domain 12511 . 



A. Effective Hamiltonians 



We first employ the Hamiltonian approach to multiple scat- 
tering |23[] . The hermitian part of the dynamics in the 
absorbing resonator is captured by a hermitian M x M- 
dimensional matrix H. We assume M ^> 1 and denote the 
level spacing in the energy range of interest as A. Loss with 
absorption rate /i is modeled by adding a non-hermitian term 
—ifit, while gain with a matching rate is obtained by inverting 
the sign of \i. We also specify an N x M-dimensional cou- 
pling matrix V between the M internal modes and TV scat- 
tering states. This matrix can be chosen to satisfy V T V = 
diag(v m ), where N finite entries v m = AM/tt describe the 
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open channels, while M — N entries v m = describe the 
closed channels. The N x iV-dimensional scattering matrix 
of the absorbing resonator is then given by 

S L (u) = 1 - 2iV(uj -ifji-H + iV T V)~ V T (41) 

Thus, the leakage from the system effectively adds an addi- 
tional non-hermitian term — iV T V to the Hamiltonian. For 
VT symmetry the scattering matrix of the amplifying res- 
onator follows as l24ll 

S r (lj) = [S^iu*)]* = l-2iV{u+ifj,-H*+iV T V)- 1 V T . 

(42) 

We also include a semitransparent interface, with scattering 
matrix (l34l) . For a closed system, the scattering quantization 
condition (l40b can then be rearranged into an eigenvalue prob- 
lem det (uj — %) = with effective Hamiltonian [12411 

*-(W r + *)- (43) 

Here Y = diag (jm) is a rea l positive semi-definite coupling 
matrix related to V T V, but with the N non-vanishing entries 
7m + Vl - T)]AM/tt 7 modified to account 

for the finite transparency T of the interface. 

The effective Hamiltonian obeys the relation VT[H] = 
cr x 7-L*cr x = %. In a PT-symmetric basis the secular equation 
det (lj — T-L) = then takes the form of a polynomial with 
real coefficients, which guarantees that the eigenvalues are ei- 
ther real or occur in complex-conjugate pairs, as required by 
equation ©. 

Analogous considerations apply to VTT' - symmetric sys- 
tems l26D • The scattering matrix of the right resonator 
then changes to 

Sr(lj) = KV*)] 1 ' = l-2iV(uj^ifi-H+iV T V)- 1 V T , 

(44) 

and the effective Hamiltonian of the closed system takes the 
form 

We then have VTT'[H] = cr x 1-0cr x = H, which again guar- 
antees the required spectral constraints. 

While Hamiltonians with these symmetries could simply be 
stipulated, the derivation of the specific manifestations (1431) 
and (l45l) reveals further constraints dictated by the coupled- 
resonator geometry. In particular, the coupling is only phys- 
ical if the matrix F is positive semidefinite. These Hamilto- 
nians bear a striking resemblance to models of mesoscopic 
superconductivity 11191 13111 . 

B. Quantum maps 

Alternatively, one can set up effective descriptions in the 
time domain B25I1 . For VT symmetry this sets out by writing 
the scattering matrices as D3I435I1 

S l (uj) = WF[exp(-iuT + n) - QF]- 1 ^, (46a) 

Sr(lj) = WF T [exp(-icur - fj) - QF T ]- 1 W T , (46b) 



(b) 
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FIG. 4: (Online version in colour.) Interpretation of the VT- 
symmetric quantum map d48l ), which translates the resonator dynam- 
ics in (a) to a stroboscopic evolution in two coupled Hilbert spaces 
(b). For the VTT 1 - symmetric map (50), F T is replaced by F. See 
also ||2l. 

where F is an M x M-dimensional unitary matrix which 
can be thought to describe the stroboscopic time evolution 
between successive scattering events off the resonator walls. 
The transfer across the interface between the two resonators 
is again described by an N x M-dimensional coupling ma- 
trix W. However, the combination P = W T W now projects 
the wave function onto the interface, such that rank P — N 
is the number of channels connecting the resonators, while 
Q = 1m — P is the complementary projector onto the res- 
onator wall (rank Q = M — N). The stated frequency depen- 
dence corresponds to stroboscopic scattering with fixed rate 
t -1 ; the variable uo thus plays the role of a quasienergy. The 
parameter /i > again determines the absorption and ampli- 
fication rate. For the passive system (/i = 0), the scattering 
matrices are unitary, which corresponds to the hermitian limit 
of the problem. 

With these specifications, and including a semitransparent 

interface parameterized by a = \J \[R + i\/T, the quantiza- 
tion condition (f40b is equivalent to the eigenvalue problem 

Tipn = A n ?/Vi, A n = ex.p(-iuj n r) (47) 
for the quantum map 

^ C ={ -UmaP ReaP + Q ) (48) 
where we symmetrized the coupling 

The 2M x 2M dimensional matrix T can be interpreted 
as a stroboscopic time evolution operator acting on 2M- 
dimensional vectors ip = Q^), where and i/jr give the 
wave amplitude in the absorbing and amplifying subsystem, 
respectively. This is illustrated in figure @] The VT sym- 
metry of the quantum map manifests itself in the relation 
T = (TxlJ 7-1 ]*^, which parallels the symmetry (f20l) for the 
scattering matrix. 

Analogously, VTT' symmetry results in the quantum map 

T =^{ e ~7 e^)^- (50) 

This obeys the symmetry T = cr^JJ 7-1 ] ^a x , which parallels 
equation ([26b . 
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The spectral properties associated with these symmetries 
are now embodied in the secular equation s (A) = det (F — 
A) =0, which exhibits the mathematical property of self- 



inversiveness 



,(1/A*) = [A- 2M S (A)]* 5 (0), 



(51) 



where s(0) = det J 7 = (detF) 2 . For each eigenvalue A n , 
we are thus guaranteed to find the eigenvalue An = [A" 1 ]* = 
exp(— zcj*t), which recovers the constraint ©. 



X. MESOSCOPIC ENERGY AND TIMES SCALES 

In mesoscopic physics, a large range of spectral, thermo- 
dynamic and transport phenomena are fully characterized by 
a few universal time and energy scales [19, 20, 37]. We now 
have all the tools to make contact to these concepts and pro- 
vide a phenomenological description of the effects of multiple 
scattering on the spectral features of the considered coupled 
resonators. These effects generally set in when a large number 
M ^> 1 of internal modes in an energy range MA is mixed by 
scattering with a characteristic time scale r = h/MA that is 
much less than the dwell time in the amplifying or absorbing 
regions, r <C tdweii- The coupling strength between these 
regions is characterized by the associated Thouless energy 
Et = h/t dwell ~ NT A, which can be small or large com- 
pared to the level spacing A, depending on the value of the 
dimensionless conductance g = NT of the interface. Non- 
hermiticity adds the new scale /i, whose interplay with the 
other scales determines the transition from real to complex 
eigenvalues. 

In the mesoscopic regime of M ^> N ^> 1 this transi- 
tion can be expected to be universal, with the features only 
depending on M/N ^> 1, M ^> 1, and T, which we con- 
sider fixed by the geometry, as well as the variable fi. The 
transition can then be investigated by combining the effective 
models set up in the previous section with random-matrix the- 
ory [0j|[22|,[38|], thus, ensembles of Hamiltonians H (usually 
composed with random Gaussian matrix elements) or time- 
evolution operators F (distributed according to a Haar mea- 
sure) which are only constrained by the symmetries of the 
problem. Closer inspection identifies two natural scenarios 
1124 12611 . described in the following two subsections, and a 
semiclassical source of corrections to random-matrix theory, 
discussed thereafter l25ll . The interplay of the various time 
and energy scales is illustrated in figure [5j 



characteristic energy scales 
f »o 



NTA=E^ 



real 

spectrum 



complex 
spectrum 



characteristic time scales 



strongly 
►amplified 

states 



MA 
— \— 



l/fi 



dwell 



onsetting qm-cl correspondence 



= X A In M 



FIG. 5: (Online version in colour.) Sketch of the energy and time 
scales governing the spectral features of VT- symmetric resonators 
with separate T' invariance (no magneto-optical effects), as her- 
miticity is broken with absorption and amplification rate \i. This uni- 
versal picture holds under the assumption of complex wave dynamics 
(M > N, NT > 1), where M = h/rA is the number of internal 
modes mixed by multiple scattering with scattering time r and N is 
the number of channels coupling the resonators, with transparency 
T. Here A is the mean level spacing, tn — h/A the Heisenberg 
time, and Et = V^dweii = NT A the Thouless energy associated 
with the dwell time tdweii in each resonator. The critical rate at which 
eigenvalues turn complex is /io — \fNA/2n. When the Ehrenfest 
time tEhr = A -1 In M (with Lyapunov exponent A) becomes com- 
parable to the dwell time tdweii quantum-to classical-correspondence 
sets in and suppresses multiple scattering, which reduces the number 
of strongly amplified states. See also 1124 12511. 



then finds 



H + T ifi 
i/i H-iT 



(52) 



which reveals the emerging T symmetry in the hermitian 
limit. The same structure arises for systems with magneto- 
optical effects which only display VTT' . For both cases, at 
/i = the system decouples into two independent sectors. For 
T = the level sequences of these sectors coincide. For very 
weak coupling (T « 1/iV e T c , thus g <C 1) the degen- 
eracy is slightly lifted, and one can apply almost-degenerate 
perturbation theory to identify the scale /i w N^/TA/2tt 
at which typical eigenvalues turn complex. As T exceeds 
1/N = T c (g > 1), however, the levels of the original se- 
quences cross, and when one then increases /i the bifurcations 
to complex eigenvalues occur between levels that were origi- 
nally non-degenerate. In this regime, a macroscopic fraction 
of complex eigenvalues appears on a coupling -independent 
scale [i « /io = \[N Aj2i\. 



A. Systems with coupling-driving level crossings 

In ordinary optics with exact T' symmetry, systems with 
VT symmetry also display VTT' symmetry. The internal 
Hamiltonian H = H* = H T is then real and symmetric, 
and the effective Hamiltonians (l43t and (1451) coincide. (Anal- 
ogously, F = F T , so that the quantum maps (1481) and (f50l) 
coincide as well.) Switching to a parity-invariant basis one 



B. Systems with coupling- driving avoided crossings 

For systems with VT symmetry but broken T' symmetry, 
the parity basis does not partially diagonalize the effective 
Hamiltonian. At fi = T = 0, one still starts from two de- 
generate level sequences, but these interact as T is increased, 
and instead of level crossings one observes level repulsion. 
In this case the crossover to a complex spectrum appears at 
ji w y/NTA 1 2ir and thus is always coupling dependent. 
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C. Growth and decay rates in the semiclassical limit 

In the semiclassical limit of large M at fixed M/N the re- 
sults above imply that the crossover to a complex spectrum 
appears at a vanishingly small rate /i <C E T . Keeping n/E T 
in this limit fixed and finite, the real phase is thus completely 
destroyed, and the focus turns to the typical decay and growth 
rates encoded in the imaginary parts Imcj n of the complex 
eigenvalues. Numerical sampling of the random-matrix en- 
sembles suggests that the distribution P(lm uj n j \i\ /jl/Et) at- 
tains a stationary limit for M > 1 l24l - [26ll . At /i > E T , 
one then finds a finite fraction of strongly amplified states 
with Imco> n & /i, thus, a large number of candidate lasing 
states. However, the fraction of these states is reduced when 
one takes dynamical effects into account 112511 . One then finds 
that the strongly amplified states are supported by the clas- 
sical repeller in the amplifying parts of the system, whose 
fractal dimension dn is more and more resolved as the phase 
space resolution h oc 1/M increases in the semiclassical limit. 
This phenomenon can be characterized by the Ehrenfest time 
tehr = A -1 InM, where A the Lyapunov exponent in the 
classical system. If ^Ehr > ^dweib the wave dynamics be- 
come quasi-deterministic, and multiple scattering is reduced. 
The fraction of strongly amplified states then follows a frac- 
tal Weyl law (a power law in h with non-integer exponent), 
a phenomenon which was previously observed for passive 
quantum systems with finite leakage through ballistic open- 
ings Mm. 



XL SUMMARY AND OUTLOOK 

In summary, the scattering approach proves useful to de- 
scribe general features of non-hermitian optical systems with 
VT and VTT' symmetry. In particular, the approach fully 
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